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Statement  of  Problem 


This  has  been  a  three  year  study  in  the  investigation  of  topics  in  the  general  area 
of  statistical  computing.  For  example,  joint  work  with  Dr.  Malcolm  Taylor  at  BRL  has 
led  to  the  SIMDAT  algorithm  for  using  a  set  of  projectile  data  to  generate  other  data 
sets  which  could  have  happened.  Also,  we  have  developed  a  simulation  based  tech¬ 
nique,  SIMEST,  to  handle  problems  not  tractable  via  the  closed  form  maximum  likeli¬ 
hood  approach.  This  technique,  whose  applicability  we  have  amply  demonstrated  on 
oncological  models  where  classical  procedures  had  stalled,  is  now  being  examined  in 
econometric  contexts.  We  have  also  developed  epidemiological  models  of  AIDS  and 
discovered  that  the  epidemic’s  sudden  appearance  is  very  likely  due  to  the  enhance¬ 
ment  factor  provided  by  bathhouses.  We  have  built  our  earlier  work  in  density  estima¬ 
tion  to  give  automated  procedures  for  selecting  bandwidths  in  the  kernel  estimation 
approach.  We  have  developed  graphical  techniques  for  demonstrating  density  esti¬ 
mates  in  dimensions  as  large  as  6.  We  have  done  extensive  work  to  build  algorithms 
whose  function  is  the  creation  of  multicolored  graphical  displays  of  density  estimates. 
We  have  done  work  in  parameter  estimation  for  time  series  problems  using  both  fre¬ 
quency  based  and  time  based  approaches.  We  have  also  done  work  for  using  in  the 
general  smoothing  of  physical  processes. 


Summary  of  Work 


July  1,  1985-December  31,  1985 

The  work  of  Chiu  has  been  concerned  with  the  problem  of  estimating  the  variable 
time  delays  of  a  signal  arriving  at  an  array  of  sensors.  A  procedure  to  estimate  the 
parameters  of  a  linear  time  delay  model  has  been  developed.  The  procedure  compares 
the  Fourier  transforms  at  different  frequencies  (thereby  taking  the  Doppler  effect  into 
consideration).  Under  regularity  conditions,  the  estimate  obtained  is  shown  to  be  con¬ 
sistent  and  asymptotically  normal.  Simulations  were  carried  out  and  the  results  were 
found  to  agree  well  with  the  theoretical  results.  The  procedure  was  applied  to  the 
records  of  the  Imperial  Valley  earthquake  of  October  15,  1979  as  recorded  by  the  El 
Centro  differential  array. 

Also,  Chiu  has  considered  the  problem  of  smoothing  noisy  images.  His  pro¬ 
cedure  preserves  some  important  characters  of  the  images,  such  as  edges  and  lines. 
The  smoothed  images  can  then  be  used  to  detect  and  locate  the  edges  of  the  images. 

Scott  has  continued  his  investigations  of  practical  and  theoretical  issues  in  non- 
parametric  density  estimation.  He  has  introduced  new  powerful  density  estimation 
algorithms  for  multivariate  estimation.  He  has  examined  the  theoretical  asymptotics  of 
several  important  density  estimators.  Algorithms  for  automatically  calibrating  these 
density  estimators  have  been  introduced  and  examined  using  asymptotic  analysis  and 
Monte  Carlo  techniques.  Practical  experiences  with  these  methods  have  been  reported 
in  several  papers  and  a  chapter  of  a  book. 

Thompson  has  been  concerned  with  the  development  of  simulation  based  algo¬ 
rithms  for  parameter  estimation  in  applied  stochastic  process  models.  Also,  he  has 
continued  the  application  of  such  models  to  oncological  situations.  He  was  one  of  four 
statisticians  selected  to  overview  recommendations  of  EPA  scientists  to  institute 
draconian  controls  on  sulfur  dioxide  emissions  from  factories  in  Ohio  and  other  states; 
this  review  found  that  such  a  move  was  contraindicated  by  the  facts;  the  proposed  con¬ 
trols  were  shelved.  An  historical  analysis  of  the  development  of  quality  control  in  the 
United  States  has  been  used  as  a  means  of  analyzing  current  difficulties  and  proposing 
solutions. 


January  1,  1986  -  June  30,  1986 

The  doctoral  dissertation  of  Husseman  demonstrates  the  use  of  variable  tile  width 
histograms  for  the  nonparametric  estimation  of  two  dimensional  densities.  It  is  shown 
how  the  procedure  can  be  used  to  obtain  robust  estimates  in  a  variety  of  situations. 

Scott  has  examined  several  classical  cross  validation  techniques  for  density  esti¬ 
mation.  An  integrated  and  highly  automated  generalization  of  earlier  techniques  is  pro¬ 
posed. 
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A  procedure  for  using  limited  data  to  generate  simulated  outcomes,  SIMDAT,  in 
a  tank  battle  is  proposed  and  demonstrated  by  Taylor  and  Thompson.  The  procedure 
represents  a  significant  improvement  over  the  bootstrap  in  many  situations. 

The  book  Cancer  Modeling  is  the  first  to  bring  modem  statistical  computer  based 
modeling  and  simulation  techniques  together  for  use  in  oncological  situations. 

The  SIMEST  algorithm  of  Thompson  et  al  enables  the  user  to  go  directly  from 
micro-axioms  to  parameter  estimation  without  attempting  to  obtain  a  “closed  form.” 


July  1,  1986  -  December  31,  1986 

The  work  of  Chiu  was  concerned  mainly  with  improved  techniques  for  parameter 
estimation  of  time  series  in  the  frequency  domain.  Additional  work  was  also  carried 
out  for  the  design  of  tests  for  the  detection  of  periodic  components  in  a  white  noise 
time  series. 

Scott  prepared  200  pages  of  handouts  for  his  shortcourse  (given  jointly  with  Ber¬ 
nard  Silverman  of  the  University  of  Bath)  in  nonparametric  density  estimation  at  the 
national  meetings  of  the  American  Statistical  Association.  He  also  examined  the  possi¬ 
bility  of  the  construction  of  expert  systems  for  the  design  of  smoothing  parameters  of 
nonparametric  density  estimators. 

Thompson  wrote  a  book  for  his  shortcourse  at  the  Army  Design  of  Experiments 
Conference  at  Monterey,  which  was  distributed  to  the  participants  in  spiral  bound 
form.  Demonstration  software  for  the  shortcourse  was  programmed  by  Kauffman  and 
Macintosh  diskettes  made  available  to  the  participants.  Thompson  also  created  an  epi¬ 
demiological  model  for  AIDS,  which  arrived  at  the  conclusion  that  the  disease,  in  the 
United  States,  was  fragile  and  dependent  on  the  presence  of  a  highly  promiscuous  sub¬ 
group  of  the  homosexual  population.  The  model  presented  predictions  for  the  course 
of  the  disease. 


January  1,  1987  -  June  30,  1987 

Kauffman  has  been  working  on  a  dissertation  (Thompson,  thesis  advisor)  on 
simulation  based  techniques  for  parameter  estimation  in  stochastic  processes. 

Chiu  has  been  working  on  procedures  to  determine  the  existence  of  signal  sinuo- 
soidal  waves  embedded  in  a  white  series.  This  problem  is  important  in  signal  process¬ 
ing  the  records  of  an  array  of  sensors  (spatial-temporal  processes)  when  one  is 
interested  in  detecting  signals  propagating  through  the  array.  He  is  also  extending  his 
one  dimensional  procedure  for  smoothing  images  to  the  two  dimensional  case. 

Scott  is  pursuing  his  work  in  graphical  based  techniques  for  the  estimation  of 
densities  in  dimensions  of  three  and  greater.  He  continues  to  work  closely  with  the 
DOD  and  NASA  in  the  analysis  of  real  world  data  sets. 

Thompson  is  working  on  density  estimation  in  the  nongraphical  mode.  He  is  also 
working  on  simulation  based  techniques  for  parameter  estimation.  His  October  Army 
shortcourse  notes  in  model  building  have  been  accepted  by  John  Wiley  and  Sons  as 


the  basis  for  a  book;  he  is  currently  at  work  on  putting  the  notes  into  book  form. 

July  1, 1987  -  December  31,  1987 

Kauffman  continues  his  work  on  a  dissertation  extending  the  SIMEST  algorithm 
of  Thompson  et  al  for  parameter  estimation  in  applied  stochastic  processes  without  the 
necessity  of  finding  closed  form  solutions. 

Chiu  has  been  extending  his  smoothing  procedure  for  two  dimensional  images. 
Also,  he  has  been  working  on  procedures  to  determine  the  existence  of  signal  sinuo- 
soidal  waves  embedded  in  a  white  series.  This  problem  is  important  in  signal  process¬ 
ing  the  records  of  an  array  of  sensors  (spatial-temporal  processes).  He  is  also  extend¬ 
ing  his  one  dimensional  procedure  for  smoothing  images  to  the  two  dimensional  case. 

Scott  is  extending  his  average  shifted  histogram  technique  for  the  estimation  of 
densities  of  higher  dimension.  He  also  is  concerned  with  cross-validation  procedures 
for  his  density  estimation  algorithms.  He  works  closely  with  the  DOD  and  NASA  in 
the  analysis  of  real  world  data  sets. 

Thompson  has  edited  a  book  Cancer  Modeling  with  Barry  Brown  in  which 
modem  data  analytical  procedures  are  developed  for  dealing  with  the  problem  of 
cancer  origination,  progression  and  control.  The  SIMEST  algorithm  is  presented  in 
that  book.  Thompson  has  written  a  paper  which  argues  that  the  current  AIDS  epidemic 
would  not  have  occurred  if  public  health  authorities  had  not  permitted  the  licensing  of 
gay  bath  houses  in  major  American  cities.  Thompson’s  1986  Army  shortcourse  notes 
in  model  building  have  been  accepted  by  John  Wiley  and  Sons  as  the  basis  for  a  book; 
it  is  hoped  that  the  final  manuscript  can  be  given  to  the  publisher  shortly. 

January  1, 1988  -  June  30,  1988 

Kauffman  is  extending  the  SIMEST  algorithm  of  Thompson  et  al  for  parameter 
estimation  in  applied  stochastic  processes  without  the  necessity  of  finding  closed  form 
solutions.  The  extension,  which  deals  with  the  inclusion  of  covariate  data,  will  be 
Kauffman’s  doctoral  dissertation. 

Chiu  continues  his  frequency  domain  approach  to  time  series  analysis.  He  has 
been  dealing  with  smoothing  (bandwidth)  considerations.  He  is  also  examining  fre¬ 
quency  domain  based  techniques  for  parameter  estimation. 

Scott  has  been  engaged  in  an  collaborative  investigation  with  Dr.  Wolfgang  Haer- 
dle  of  the  University  of  Bonn.  This  has  involved  visits  by  Scott  to  Bonn  and  by  Haer- 
dle  to  Rice.  The  work  involved  has  included  topics  in  robust  regression  and  in  mul¬ 
tivariate  nonparametric  density  estimation.  The  utilization  of  these  techniques  in 
mathematical  economics  has  been  a  major  thrust  of  their  work. 

Thompson  has  been  examining  multiresponse  variable  versions  of  the  SIMEST 
algorithms.  The  copy  editing  step  of  the  publication  of  Thompson’s  book  Empirical 
Model  Building  has  been  completed  by  John  Wiley  &  Sons.  The  target  date  for  publi¬ 
cation  is  the  end  of  1988.  Thompson  has  consulted  with  scientists  at  Fort  Ord  on  the 


subject  of  data  based  war  gaming. 
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